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Over the past few years one of us (Murthy) in collaboration with R. Shankar has developed an 
extended Hamiltonian formalism capable of describing the ground state and low energy excitations 
in the fractional quantum Hall regime. The Hamiltonian, expressed in terms of Composite Fermion 
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operators, incorporates all the nonperturbative features of the fractional Hall regime, so that con- 
ventional many-body approximations such as Hartree-Fock and time-dependent Hartree-Fock are 
applicable. We apply this formalism to develop a microscopic theory of the collective edge modes in 
fractional quantum Hall regime. We present the results for edge mode dispersions at principal filling 
factors v — 1/3, 1/5 and v — 2/5 for systems with unreconstructed edges. The primary advantage 
of the method is that one works in the thermodynamic limit right from the beginning, thus avoiding 
the finite-size effects which ultimately limit exact diagonalization studies. 
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I. INTRODUCTION 



The bulk fractional quantum Hall effect is by now well-understood as a consequence of an interaction-driven 
incompressible stat o 1 ' 2 ' 3 of a two-dimensional electron gas (2DEG) in a strong perpendicular magnetic field. It was 
also realized long ago^ that the edges of the bulk system play a fundamental role in transport at low frequencies, since 
V ! the bulk is incompressible and the gapless excitations are available only at the edges. An effective description of edge 
excitations for incompressible fractions was first proposed by Weni£ He developed a hydrodynamic theory of a sharp 
edge, which is a realization of the one-dimensional Chiral Luttinger liquid. 

Recently, several high-precision experiments have probed the low-frequency dynamics of the edge in the fractional 
quantum Hall (FQH) regime over a range of filling factors&i Some features of the results, such as the observation of 
a power law in the I-V curve which depends on the filling factor, and the magnitude of the power-law exponent^ are 
unexpected from the point of view of hydrodynamic theories of the edgeA^ Many explanations have been proposed 
for these discrepancies)2ii2iiiii2ii&ii±£ii& including the effects of long-range interactions and edge reconstructions. 
Motivated by these experiments, Wan et alM have recently explored the effects of edge reconstructions on the 
behavior of the collective modes using numerical exact diagonalization. While important qualitatively new results^ 
£f) . have been obtained by this research, these studies do suffer from finite-size effects. 

We want to investigate the physics of FQH edges using the extended Hamiltonian theory developed by Shankar 
and Murthysii This approach has the advantages that it starts with a microscopic Hamiltonian, permits analytic 
calculations in the thermodynamic limit and also retains the lowest Landau level limit, ft thus provides a complemen- 
tary approach to that of exact diagonalization studies, fn a previous papeni^ we applied the extended Hamiltonian 
theory to study the structure of FQH edges within the Hartree-Fock (HF) approximation, fn this paper, we further 
explore the edge states beyond mean-field level. We use the time-dependent Hartree-Fock (TDHF) approximation 
which systematically incorporates particle-hole states because it is a conserving approximation^ i.e. it respects the 
c*j ' symmetries of the problem. 

The structure of our paper is as follows: In the next section, we review the essentials of the extended Hamiltonian 
theory and mean-field theory of FQH edges. The time-dependent Hartree-Fock approximation is described in Sec- 
tion ^Q] Our numerical results for v — 1/3, 1/5 and v = 2/5 are presented in Section ITVI The results for v = 2/5 are, 
to the best of our knowledge, the first in the literature. The last section contains conclusions and caveats. Throughout 
this paper we restrict ourselves to unreconstructed edges due to the computational limitations of our own. 



II. EXTENDED HAMILTONIAN THEORY AND HARTREE-FOCK APPROXIMATION 



For a two-dimensional electron gas in a strong perpendicular magnetic field, the kinetic energy of electrons is 
quenched and the single-particle states are organized into Landau levels separated by the cyclotron energy. For 
electrons in the lowest Landau level, neglecting the zero point energy, the interaction Hamiltonian is given by 

H * = 2l E V ^ : PMpei-Q) ■■ (1) 
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where V ee (q) is the repulsive electron-electron interaction. This interaction energy scale is typically several tens of 
Kelvins, of the same order as the cyclotron energy (about 100K). Nevertheless, in the integer quantum Hall regime, the 
interaction is usually treated as a perturbation on the ground-state Slater-determinant of v Landau levels completely 
filled. For a fractional filling factor, this ground state is degenerate and the perturbative approach fails i The key 
insight of Jain^was that if one thinks in terms of Composite Fermions (CFs), which are envisioned as electrons attached 
to an even number 2s of flux quanta, one can use the independent particle picture for the Composite Fermions. In this 
approach, a system at filling factor v = pj (2ps + 1) has CFs which see an effective field B* = B/(2ps + l), which is just 
right for them to fill p CF Landau levels (CF-LLs). Based on this insight, Jain constructed excellent wavefunctions 
for ground state and low-lying excited states for the principal filling factors. However in order to calculate response 
functions or finite-temperature properties, one needs a dynamical theory. Initial attempts focused on bosonioSS or 
fermionioSi variants of the Chern-Simons theory^ culminating 23 in a Chcrn-Simons theory for v = 1/2, which is a 
Fermi-liquid-like stated These theories, while satisfactory in many respects, still did not incorporate some of the 
nonperturbative aspects of the CF, such as the facts that in the lowest-Landau-level limit the effective mass of the 
CF is determined by the interactions alone or that its charge is fractional. 

To handle these problems, Shankar and Murthy developed the extended Hamiltonian theoryjii a detailed account of 
which can be found in a recent review^ A key ingredient of this theory is the introduction of a new set of coordinates 
for pseudovortices in addition to those of electrons. For a FQH liquid at filling factor v — p/(2ps + 1) where p and s 
are integers, each electron couples with 2s pseudovortices to form a composite fermion (CF) in the following way 

I 2 

Re = r--—zx\l, (2) 
1 + c 

I 2 

R « = f + Ti : \ 2 x n ' ( 3 ) 

c(l + c) 

where c 2 = 2vs, I — \Jh/eB and r and n are the position and velocity operators of the composite fermion. The 
electron guiding center R e and the pseudovortex guiding center R v satisfy the algebra 

[R ea ,Rv0] = 0, (4) 
[Real Rep] = —^^ap, 

r i - /2 

[Roa, Rvp\ — +i-^£al3- 

Thus, the electron guiding-center coordinates satisfy the magnetic algebra with charge — e < 0, whereas the pseu- 
dovortex guiding-center algebra represents an object with charge +ec 2 . The CFs thus have a magnetic algebra charge 
of e* = — e(l — c 2 ), showing that they are subject to a reduced magnetic field B* — Bj (2ps + 1), just right to fill the 
first p CF-Landau levels, exactly as in Jain's pictured The Hartree-Fock state of CFs with the first p Landau levels 
filled provides a nondegenerate starting point for analytical calculations. 

To calculate the matrix elements of these operators in the CF basis, we use the single-particle states of the CFs in 
the reduced effective field B*. In the Landau gauge, a single-particle state \nX) is characterized by the CF-Landau 
level index n and the guiding-center coordinate X. In the real-space representation, the (unnormalized) single-particle 
wavefunction is 

(r\nX) = e lXy ' 1 ' 2 e~^~ x)2/21 '* H n [(x - X)/l*} , (5) 



where H n (x) are the Hcrmite polynomials and /* = ^Jh/eB* is the magnetic length in the reduced field B* seen by 
the CFs. In the following, we use I* = 1 as the unit of length. Using this basis, it is straightforward to express the 
electron density p e and the pseudovortex density p v operators in second-quantized notation 

Pe{q) = ^ e ~ lq * :X dl lX _ qy/2 d n2 x+ qy /2Pn 1 n 2 {q)i (6) 
rnX 

Pv{q) = J^^^^niX-qy^^X+qy^X^nM, (7) 



jX 



where d n x($ n x) destroys (creates) a CF in the state \nX), and p ni n 2 (Xn 1 n 2 ) are the plane-wave matrix elements for 
electron guiding center R e (pseudovortex guiding center R v )w^ We have now expressed the electron density operators 
and hence the microscopic Hamiltonian for the fractional quantum Hall systems in terms of CF operators, d n x and 
d* nX . Therefore the original problem of interacting electrons in the fractional quantum Hall regime has been readily 
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transformed into a problem of CFs in the integer quantum Hall regime where various many-body approximations are 
applicable. But, we must remember that there is no free lunch; the price to pay is that the transformed problem is 
subject to special constraints as a result of the introduction of pseudovortex coordinates (section HTT|) . 

We next consider the positive neutralizing background charge Pb(x) which produces a confining potential near the 
edge at x = 0. The corresponding Hamiltonian is 

H b = ^V eb (-q)p e (<?) (8) 

Q 

where 

/>oo 

V eb (x) = dx'V eb (x-x')p b (x'), (9) 



and V eb (x) is the attractive electron-background interaction. We choose a background charge density which vanishes 
linearly over width W near the edge: p b (x) = for x < —W/2, p b (x) = — ^ [l + |p] for \x\ < W/2, and p b {x) = p a 
for x > W/2 where po is the CF density in the bulkii&Si We also assume that the background charge resides in the 
same plane as the 2DEG, V eb (q) — —V ee (q). 

In the bulk of a FQH liquid at filling factor v = p/(2ps + 1), p CF Landau levels are completely filled. Due to 
the presence of the confining potential close to the edge, different CF-LLs mix together. Since the system is still 
translationally invariant along the y-direction, the Landau gauge quantum number X is a good quantum number, 
and only CF-LLs with the same value of X mix. We introduce the mixing matrices a(X) as function of X 

d n x = a nm (X)a mX (10) 
m 

where a mX denote the annihilation operators for the CF quasiparticles which result from mixing of different CF-LLs. 
The total Hamiltonian, H e + H b , represented in terms of new operators a mX and (J mX , is decoupled using standard 
HF approximation^ and the mean-field solution is chosen to be 

(al x a nX ')=S mn S xx ,N^(X) (11) 

where N^(X) is the CF-occupation number at guiding center X and CF Landau level m. Note that with this 
ansatz for the mean-field solution we are explicitly forbidding any structure in the y-direction. Upon decoupling, the 
mean-field Hamiltonian becomes 

h hf = H mn (X)al x a nX (12) 

mnX 

We determine the mixing matrices a's so that the mean-field Hamiltonian H HF is diagonal. The process is normally 
carried out iteratively as in our previous work^ The quantities of interest in mean-field study are the energy bands 
E m {X) = H mm (X), the occupation functions N^(X), and the mixing matrices a(X). 



III. TIME-DEPENDENT HARTREE-FOCK APPROXIMATION 



A. Time-dependent Hartree-Fock Hamiltonian 



As we mentioned earlier, special care is required in treating the transformed Hamiltonian. In particular, the 
original electronic Hamiltonian knows nothing about the pseudovortex density p v and thus is invariant with respect 
to our choice of it. The transformed problem has extra degrees of freedom, the pseudovortex coordinates. Exact 
solutions to this problem will automatically be degenerate along the pseudovortex manifold, and will thus have a 
gauge symmetry. However, the mean-field ground state does not have the required gauge invariance. The situation 
is similar in the theory of superconductivity, where the Bardeen-Cooper-Schrieffer (BCS) ground state is not gauge 
invariant. Nonetheless gauge- invariant dynamical response functions can still be obtained from the mean-field ground 
state by a so-called conserving approximationii^ TDHF approximation turns out to be the simplest of the class 
of conserving approximations, and is standard in the context of superconductivity when dealing with the gauge 
symmetry. 1 - Read was the first to apply the TDHF approximation to study the problem of v — 1 bosons, 27 which 
has important similarities to the FQH system at filling factor v — 1/2.24 It was subsequently used by Murthy 2 ^ to 
find the collective modes for the incompressible fractions using the extended Hamiltonian theory. 
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In this Section, we will exploit the TDHF approximation to obtain the dispersion relations of the edge excitations. 
The way the conserving nature of the TDHF approximation manifests itself is that the spectrum of the TDHF 
Hamiltonian necessarily contains zero eigenvalues, each corresponding to an unphysical gauge degree of freedom, while 
the rest of the spectrum is physical^ It is important that the unphysical zero-modes are identified and excluded and 
we will discuss the techniques we use for the same in the next subsection. 

It follows from Eq.Q) that [pe(q), Pv(<f)] = 0. The total Hamiltonian, H e + Hb, therefore commutes with the 
pseudovortex density p v . Let us consider the magnetoexciton operator 

O mim2 {X,q y ) = a l ni ^X-q y /2 a m2,X+q v /2 (13) 

where a m x is the rotated annihilation operator related to the CF-operator d as in Eq. (|10|l . We have explicitly 
made use of the fact that there is translational invariance in the y-direction and therefore q y is a good quantum 
number. Henceforth we shall suppress the index q y and simply write O mim2 (X). The equation of motion for the 
magnetoexciton operator is 

-i^-O mim2 (X) = [H, O mim2 (X)} (14) 
Using the standard HF decoupling to evaluate the right-hand-side we get 

[H, O mim2 (X )} -> [E mi (X Q - q y /2) - E m2 (X + q y /2)]O mim2 (X ) + [N^ i Xa _ qy/2 - N£ 2)Xo+qy / 2 ] 
x [h^k'^k-jiXo - X,-q y ) - h klk2k > ik > 2 (X — X ,q y )] 

Xvkik[ 

x 4 ini - q y /2)a k2n2 (X + q y /2)al, m2 (X + q y /2)a k , rni (X - q y /2)O nin2 (X) (15) 

where 

h klk2k ' lk ' 2 (X,q y ) = J ^V ee (q)p klk2 (q)p k ' ik ' 2 (-q)e^ x (16) 

and (i = (mi, 1712), v — (^1,^2) e ^c. We define a vector ^!^(X) corresponding to the operator 

O^^^WO^I). (17) 

In the space spanned by , the TDHF Hamiltonian reads 

H.{p,iy]X 0l X;q y ) = 8 mini 8 m2n2 [E mi (X - q y /2) - E m2 (X + q y /2)] + [N* u x -q y /2 ~ N m 2 ,x +q y /2] 
x ^2 [hkik' 2 k[k2 (Xq — X, — q y ) — h klk2k ' ik ' 2 (X — Xq, q y )] 

x 4 in A X -%/ 2 ) a k2n 2 {X + qy/2)al lm2 {X + q y /2)a k2mi (X -q y /2) (18) 

The TDHF Hamiltonian is the key quantity containing information about the collective modes in the system. Upon 
diagonalization, we obtain the energy spectrum of excitations, and the left (ip L ) and right (tp R ) eigenvectors. Direct 
inspection of Eas. l|ltjfl and (|18|) shows that H changes sign, after a suitable reshuffling of indices, as q y does. That 
means the spectrum at a q y contains information about the spectrum at —q y as well. The left- and right-eigenvectors 
are identical (modulo reshuffled indices) for both ±q y . We are interested in the gapless excitations across the edge. 
These modes are in the regime of low energy and momentum, and our task at hand is to identify these modes from 
the energy spectrum. 

B. Identification of the Physical Modes 

Due to the gauge symmetry of the original Hamiltonian the TDHF Hamiltonian 7i necessarily has zero eigenvalues, 
even for a nonuniform system. In the exact solution of the original Hamiltonian each energy level is infinitely degenerate 
since it costs no energy to rotate among gauge-connected eigenstates. The zero modes of TL are thus unphysical and 
we have to discard them from its spectrum. We stress that the spurious modes have precisely zero energy only when 
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the full set of infinite CF-LLs are taken into account. Obviously, in computations, though one tries to include as 
many CF-LLs as possible, that number cannot be infinite. Therefore, the technical problem of how to separate the 
unphysical modes from the low-lying physical ones arises. There are two different ways to achieve this. The first one 
is to monitor the behavior of the energy of a mode. As more CF-LLs are added, the energy of an unphysical mode 
gets closer to zero while that of a physical mode becomes stable. In the bulk case the unphysical modes were singled 
out successfully using this procedure^ The second is to resort to the electron density correlation function. Since the 
electron density operator is gauge- invariant, unphysical modes should decouple from its correlation functions as more 
CF-LLs are added. 

In the edge case, the first method is not practical because the energy spectrum of Ji contains every mode at a given 
q y including many unphysical modes which have high q x and therefore require a very large number of CF-LLs for 
their detection. We shall instead use the second method. We define the electron density correlation function 

S(X ,XQ',q y ,w) = -i(p e (X , q v ,uj)p e (Xo, -q v , -u)) (19) 

as the ground-state expectation of a four-fermion operator. A detailed but straightforward calculation gives rise to 
the following expression for the density correlator in terms of known quantities, 



S(X ,X^,q v ,uj) = 

ct.ki 



pX 



(X - X, g y )<g (X, q v )al iki (X - q y /2)a m2k2 (X + q y /2) 



x E^«*( X o " X', -q y )^£(X', ~q y K ik3 (X' + q y /2)a n2ki (X' - q y /2) 

,vX' 

x [N[ 3 (X' + q y /2) - Nl{X' - q y /2)} (co-E a + ir, sign{E a ))- 1 
where p M (X, q y ) is the ID Fourier transform of the density matrix 

P*(X,Qv)= I ^e^ X p,(q x ,q y ) 



(20) 



(21) 



As handy as it looks, Eq. I|2(J|I still has a minor technical problem. It describes how strongly a mode couples with 
the electron density; in other words, while it does discriminate between physical and unphysical modes, it does not 
discriminate between the bulk and edge modes. We note that physical edge modes at low momentum and energy, the 
regime of interest, should become stable as more and more CF-LLs are added. In a further refinement we exclude the 
high-momentum modes by the following filtering procedure. Let us consider, instead of the pure electron density p e 
used in Ea. (|20Jl . a filtered operator concentrated near the edge 



pUx ,q y )=T / e-^- x ^f p,(X,q y ) 



x 



(22) 



that cuts off all fluctuations across the edge with length-scales longer than If. The choice of Gaussian filtering in 
Ea. (|22|l is convenient for calculations. The results we present in the following section are based on the choice If = I* , 
which seems to work well. The filtered density correlator, S^(Xo, Xq, q y , to) — —i(p((Xo,q y ,u})p((XQ,—q y ,—uj)) 
satisfies exactly the same formula (|2(JH except that the electronic density matrix elements p^(X,q y ) are replaced by 
the filtered density ones pUX^qy), defined as 



J^e-^V^ x p,(q x ,q y ) 



(23) 



IV. RESULTS 



In the following calculations, we vary the confining potential at the edge by changing the width W over which the 
background charge vanishes linearly from its bulk value to zero. We consider a short-ranged Gaussian interaction 

V ee (q) = V e-« 2 x2 > 2 = -V eb (q) (24) 
with range A = I*, and a Thomas- Fermi interaction 

V ee q) = = ~ V *W ( 25 ) 
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with qrFl* = 0.4, which becomes truly long ranged as qtf — ► 0. The Thomas-Fermi case is considerably more 
difficult for the HF part of the calculation. As qrF decreases, it takes more iterations for the HF to converge and the 
presence of the edge makes itself felt deeper inside the bulk. We do not find any qualitative difference in the edge 
collective modes between the two cases. This being the case, we concentrate on the short-range Gaussian interactions 
for v = 1/5 and v — 2/5. 

Deep in the bulk, the quantum Hall liquid is uniform and the structure of CF-LLs and excitations in TDHF is 
well-understoodi2£ For numerical purposes we focus on the region around the edge, taking 8^* in the bulk side and 21* 
in the empty side. The lattice spacing is Aa; = 0.11*, resulting in 100 sites. Typically, 8 to 10 CF-LLs are considered. 
The TDHF Hamiltonian size is of the order of a few thousand, making the diagonalization feasible. Notice that Ax sets 
the lower bound on resolution for q y (recall that the magnetoexciton contains CF creation and destruction operators 
separated by distance q y .) We believe that our choice of Ax = 0.11* is appropriate as the natural length-scale in the 
problem is I*. Reducing Ax, while allowing us to explore smaller values of q y , would require more computational 
resources. 




FIG. 1: Typical mean-field energy bands (left) and electron density profile (right) at the edge for v = 1/3. We have used the 
Gaussian interaction, W = and 10 CF-LLs. The energy is measured in units of e 2 /el* and length is measured in units of 1* 
(Note that our interaction is not normalized). Only one CF-LL is occupied in the bulk, as signified by the bulk electron density 
p e = 1 in the right panel; the dotted line in the left panel denotes the chemical potential. 



First we present results for the v = 1/3 case, which corresponds to s = 1 = p or c = y/2/3. In this case, the ground 
state in the bulk has only the lowest CF-LL filled. We solve the HF mean-field equations to obtain the energy bands, 
the mixing matrices, and the CF occupation profile. Figure^shows typical mean-field energy bands (left) and electron 
density profile (right) near the edge for a Gaussian interaction, with W = and 10 CF-LLs. From the mean-field 
results the TDHF Hamiltonian 7i, Ea. (|18|) . is created and diagonalized. A direct inspection of 7i reveals that only 
the excitations from the lowest CF-LL to higher ones are relevant to the edge modes. We thus need to keep only 
elements in H that have mi = or m 2 = 0. The size of the TDHF Hamiltonian with 10 CF-LLs is 100 x 19 = 1, 900 
rows (and columns). While diagonalizing a matrix this size is easy, the most time-consuming part of the computation 
is the evaluation of the matrix TL due to the quadruple sum in Eq. (|18|) . Although the TDHF Hamiltonian is not 
symmetric, remarkably, our diagonalization results find that it only possesses real eigenvalues, meaning every mode 
is long-lived as normally expected. Upon diagonalization, the set of left- and right- eigenvectors are used to compute 
the filtered density correlator. 

Figure |21 shows the imaginary part of the local filtered correlator S' (0, 0; q y = 0.3, ui). It is evident from the picture 
that while there are several peaks in Im at 8 CF-LLs, many of them at low energy disappear at 10 CF-LLs. These 
peaks represent spurious modes which decouple from the electron density in the limit of infinite number of CF-LLs. 
The arrows indicate peaks corresponding to the physical modes which are stable as more CF-LLs are added. The peak 
at positive energy is the edge mode; it changes little in position and weight as the number of CF-LLs is increased. 
The gradual disappearance of unphysical modes and the stability of the physical modes with respect to increasing 
number of CF-LLs are strong indications that our results, although obtained using a finite number of CF-LLs, are 
reliable. 

By carrying out this procedure for different values of q y we obtain the dispersion relation for the v = 1/3 edge 
mode, which is shown in Fig. [21 The TDHF approximation, which takes into account the collective particle-hole pair 
states, reduces the excitation energy below its mean-field value and is shown with a dotted line. As expected, the 
edge mode is gapless, and at small momenta the dispersion is linear, u> c (q) — voq y , where vq is the cyclotron velocity 
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FIG. 2: Spectral function of the local electron density correlator, Im S*(0,0;q y = 0.3, ui) for a v — 1/3 system with Gaussian 
interaction and W = 0. The dotted (solid) lines show results with 8 (10) CF-LLs. The arrows indicate poles of physical edge 
modes. Notice that the spurious modes lying between the arrows are strongly suppressed as the number of CF-LLs is increased 
from 8 to 10, whereas the physical edge mode with positive energy remains essentially unchanged in both position and weight. 
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FIG. 3: Dispersion (left) and inverse velocity (right) of v = 1/3 edge mode for a system with Gaussian interaction, W = 0, 
and 10 CF-LLs. In the left panel, the dotted line is the naive mean-field result whereas the solid line is the dispersion obtained 
from TDHF approximation. The right panel shows the momentum-dependent inverse velocity l/v(q) which fits the heuristic 
form 12611 . The suppression of velocity at large q is due to the decreasing electron density near the edge. 



of the CF along the edge. For larger momenta, we assume a momentum-dependent velocity v(q). To see heuristically 
what the natural form of this velocity might be, note that as q y increases we are making excitations deeper into the 
empty region. In this region the electron density is smaller than that in the bulk and thus the effective field B* seen 
by the CFs must be larger. This means the velocity of excitations must be smaller. Assuming a linear variation of 
the electron density and thus the effective field near the edge (Fig. ^| we arrive at the following form, 

' ' ly (26) 



V{q) v ^ 

where vq = v(q = 0) and a is a constant. The right panel in Fig. [3] shows the inverse CF velocity corresponding to 
the dispersion shown in the left panel. Given that there is only one free parameter a, the data seem to fit our ansatz 
(|26J) very well. 

Now we examine how the collective mode dispersion changes as the width of the confining potential W is increased. 
We remind the reader that as W increases, the v = 1/3 edge becomes susceptible to reconstruction irrespective 
of the range of electron-electron interactionpi&Sii Therefore we expect a substantial softening of the collective mode 
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dispersion as a precursor to the edge reconstruction. Figure 0] shows the dispersion (left) and the inverse CF velocity 
(right) for W — 1.2Z*; indeed, as expected, the dispersion for W — 1.21* is softer than that for W = (Fig. EJ. 



v=1/3, Gaussian Interaction 
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FIG. 4: Dispersion (left) and inverse CF velocity (right) of the v = 1/3 edge mode with Gaussian interaction, 10 CF-LLs and 
W = 1.21*. Note that the edge mode softens as W increases (Fig. 0. 



To illustrate the effect of the range of interaction on the edge modes, we also study the v = 1/3 case with Thomas- 
Fermi interaction. In this case, it is harder to filter out the spurious modes than in the case with short-ranged Gaussian 
interactions; however, our scheme still works reasonably well. Figure [Sj shows the edge-mode dispersion (left) and 
inverse CF velocity (right) with W = and 10 CF-LLs. These results are qualitatively similar to those obtained using 
Gaussian interaction. The overall change in the edge-mode energy-scale is due to our use of un-normalized Gaussian 
interaction and has no fundamental significance. 



v=1/3, Thomas-Fermi Interaction 
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FIG. 5: Dispersion (left) and inverse CF velocity (right) for v = 1/3 state with Thomas- Fermi interaction, W — 0, and 10 
CF-LLs. Note that the dispersion is similar to that in the case of Gaussian interactions and our ansatz 12611 for inverse velocity 
works for long-ranged interactions as well. 



Next we consider filling factor v — 1/5 which corresponds to CFs (consisting of electrons with 2s = 4 flux quanta) 
filling the lowest CF Landau level {p — 1 and c = ^/4/5). Since there is no qualitative difference between collective 
edge modes obtained from Gaussian or Thomas-Fermi interactions, we concentrate on the Gaussian interaction which 
is more amenable to analytical calculations. Figure shows the dispersion of the edge mode (left) and the inverse CF 
velocity (right), with W = and 10 CF-LLs. We note that the edge- mode energy scale is similar to that observed in 
the v = 1/3 case. As expected, we find one edge mode and that mode softens with increasing W. We also find that 
the inverse CF velocity fits our heuristic form i|26|l and the free parameter a has value similar to that in the v = 1/3 
case. 



v=1/5, Gaussian Interaction 



v=1/5, Gaussian Interaction 



■ W=u 



1/v=57.6+51.1 q I* 



0.000 * ' 1 ' 1 ' 1 ' 1 1 ' 1 1 1 1 1 ' 1 

0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

q y l* q y l* 

FIG. 6: Dispersion (left) and inverse CF velocity (right) for the v = 1/5 edge collective mode with Gaussian interaction, W = 
and 10 CF-LLs. The fit of the inverse CF velocity (right) to our heuristic form is still reasonable. The large scatter at small 
momenta arises from uncertainties in locating the physical edge mode at small energies. 



Let us now turn to the case of a fully spin-polarized v = 2/5 state. Here we present, to the best of our knowledge, 
the first microscopic treatment of collective edge modes for v = 2/5. In this case, the CFs (consisting of electrons 
with 2s = 2 flux quanta) fill p — 2 CF Landau levels in the bulk, which leads to more involved numerical calculations. 
One now needs a TDHF matrix with twice the dimension as in the case of v = 1/3 or v — 1/5, since there are two 
filled CF-LLs in which the hole can be created. Consequently, we expect two gapless modes at the edge propagating 
in the same direction. Due to computational constraints, we restrict ourselves to Gaussian interaction and 8 CF-LLs. 
Figure [7| shows the spectral function of the local density correlator, Im S^(0, 0; q y , oj) for q y l* = 0.4 and q y l* = 0.5. 
For each value of q y , two peaks marked by arrows represent the physical edge modes. It is interesting to note that 
contrary to the naive speculation, the edge mode with smaller energy has a greater spectral weight and therefore 
higher charge. Right panel in Fig. [S] shows dispersions for both edge modes, which have approximately the same 
velocity; the left panel shows that the inverse CF velocities can be fit to the heuristic form obtained earlier. 




FIG. 7: Spectral function of the local electron density correlator for v = 2/5 with Gaussian interaction, W — and 8 CF-LLs. 
The positions of the physical edge modes, defined by their stability with increasing CF-LLs and their positive energies, are 
marked with arrows for q v l* — 0.4 (solid) and q y l* = 0.5 (dotted). 
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v=2/5, Gaussian Interaction 

• 1st edge-mode 

* 2nd edge-mode 
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FIG. 8: Dispersions (left) and inverse CF velocities of edge modes for a v = 2/5 system with Gaussian interaction, W = and 
8 CF-LLs. The two modes have similar velocities and their inverse CF velocities fit the ansatz I26H reasonably well. 



V. DISCUSSION 



In this paper we have presented a microscopic approach to calculating collective edge excitations in the fractional 
quantum Hall regime which is complementary to the exact diagonalization approach. This approach is based on the 
extended Hamiltonian theory 25 and permits the use of many-body approximations which are generally valid in the 
integer quantum Hall regime. We have presented results for collective modes, and their dependence on the nature 
of electron-electron interaction and the width of the background confining potential at various filling factors. For 
an unreconstructed edge, at filling factor v = 1/3 and v — 1/5, we found a single linearly dispersing edge mode 
which softens with the increasing width W of the confining potential. In contrast, we found two linearly dispersing 
modes for fully spin-polarized v = 2/5 state. It is somewhat surprising that the two edge modes have nearly the 
same energy, and that the mode with the lower energy has a greater overlap with the charge operator. The curvature 
of the collective mode dispersion is somewhat unexpected from hydrodynamic theories. We show that it can be 
understood satisfactorily by assuming a momentum-dependent collective mode velocity, Ea. l|26|l . which reflects the 
effect of decreasing electron density and the subsequent increase in the effective field B* seen by the CFs. Our ansatz 
implies that the magnetoexciton frequency Lo c {q) = v{q)q y depends only on the combination voq y . Incidentally, we 
find that the free parameter a is roughly equal for v = 1/3 dispersion and both branches of v — 2/5 dispersion; 
therefore we can collapse all three dispersions (which have s = 1) onto a single curve by scaling with vq (Fig- El- 
Our results for edge-mode dispersions did not show the roton minimum near q y l* ~ 1 seen in exact diagonalization 
studies, 14 presumably because the system is not close to edge reconstruction^ 

Now let us mention some caveats. While our approach does not suffer from the computational limitations of exact 
diagonalization, we have computational limitations of our own. We work, in principle, in the thermodynamic limit. In 
practice, we are forced to truncate the filled states in the bulk and the empty states near the edge, to use a discretized 
Landau gauge index X, and to use a finite number of CF-LLs in the calculation. There are two limitations which we 
face: At the mean-field stage, we need to find the HF ground state as accurately as possible. Each step of the iterative 
HF calculation takes longer when these numbers get larger. In the TDHF stage, the main limitation is imposed by 
the time required to construct the TDHF matrix which involves a quadruple sum and increases as ./V 4 where N is 
a measure of the number of single-particle levels kept. The diagonalization of the TDHF Hamiltonian takes a small 
fraction of the total computation time. We originally intended to study the softening of the collective edge-modes 
prior to edge reconstruction, and the evolution of edge-modes thereafter, as well. However, we were not successful 
in identifying the physical modes unambiguously at larger values of W . Presumably, with better computational 
resources, if we increase the number of CF-LLs used, this task will be feasible. 

What our approach does is to bring within computational reach the techniques which have been successful in the 
study of edge modes of integer quantum Hall systems^SS It does have its limitations; despite these limitations, we 
believe it can be usefully employed for cases which cannot be studied microscopically by any other method. For 
example, if one wants to calculate the temperature dependence of edge mode dispersion, one requires the knowledge 
of all excited states in an exact diagonalization approach. This is computationally prohibitive compared to finding 
the ground state and a few excited states. The state of the art exact diagonalization studies can treat at most 7-9 



11 



Scaling of ro (q) with v q I* 



0.04 | . 1 . 1 . r- 



K r A 

3 M 



• v=1/3 (single branch) 

* v=2/5 (lower branch) 

♦ v=2/5 (upper branch) 



0.00 0.02 0.04 0.06 0.08 



FIG. 9: Scaling of the single-branch of v = 1/3 dispersion and both branches of the v = 2/5 dispersion with vo. Since both 
have s = 1, and similar values of the free parameter a which appears in Eg . 1261 . the data collapse on a single curve. 



particles if all the states are keptii* We expect the collective edge modes to be quite sensitive to temperature and 
finite size effects, since they are gapless. This may be the reason that although exact diagonalization studies with 
all states^ reproduce the qualitative features of microwave absorption experiments^ the energy scales predicted are 
considerably higher than those in the experiments. Another important piece of physics, which our approach can easily 
handle, is the role of disorder. In any real sample there is bound to be disorder. If the dominant disorder is due 
to quenched fluctuations in the remote dopant layer, the system breaks up into incompressible strips separated by 
compressible regions^ These strips have been seen in imaging experiments It is straightforward to include the 
disorder potential perturbatively in our formulation, both in the mean- field and the TDHF approximations. The HF 
matrix is somewhat more complicated due to possible spatial structure along the edge. However, using the clean 
single-particle and collective modes of the incompressible strip as a basis, one can construct a microscopic theory of 
collective edge modes with disorder. This will provide us with a microscopic approach for examining the parameters 
which enter hydrodynamic theories to this problem. 8 

One can also address novel edge-states with our approach, such as an edge between v = 1/3 and v = 2/5 regions 
which have CFs with the same number of flux quanta or an edge between v = 1/3 and v = 2/9 regions, which have 
CFs with different numbers of flux quanta attached. It would also be very interesting to study the edge dynamics of 
the Fermi-liquid-like state at v = l/2m£ 
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